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Abstract 

Qh! We explore the feasibihty of directly extracting the large-x valence d/u ratio 

^ I through a measurement of pions in the current fragmentation region of semi- 

inclusive deep-inelastic scattering from protons. 

><" 
H , 

5^ I The d/u ratio contains important information about the flavor structure of the proton. 

In particular, its asymptotic behavior at large x can tell us which mechanism is responsible 
for the breaking of SU(2)spinX SU(2)flavor symmetry. Given that there are firm predictions 
for this behavior, it constitutes a serious test of perturbative QCD 

So far, a direct measurement of the d/u ratio has been rather difficult, mainly because 
the cross sections decrease rapidly in the extreme kinematics near x ~ 1. Previous analyses 
have used inclusive deep-inelastic scattering (DIS) data on proton and deuteron targets 
to obtain the neutron to proton structure function ratio, from which d/u can be extracted 
at large x according to: 

Fl_l + Adlu 

" 4 + d/u' ^ ^ 

However, the neutron structure function in Eq.(|^) is obtained from data on the deuteron, 
which suffer from the fact that nuclear effects, even in the deuteron, become quite significant 
at large x In particular, whether one corrects for Fermi motion only, or in addition for 
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binding and nucleon off-shell effects, the extracted neutron structure functions for x ^ 0.7 
can differ rather dramatically 0. 

The question is therefore how to avoid the problem of the uncertainty in the extraction 
procedure introduced by nuclear effects. A number of suggestions how to overcome this 
problem have been raised in the literature. Perhaps the simplest idea is to use neutrino 
and antineutrino beams to scatter on proton targets, thereby selecting the d and u quark 
flavors, respectively. Although such data exist []^, their statistical accuracy is unfortunately 
not very high. Furthermore, the range of x covered {x ^ 0.6) is too small to constrain the 
asymptotic behavior. Another recent suggestion was to use W^-boson production in pp and 
pp collisions [^, measuring the decay muon asymmetry at large lepton rapidities. While 
simple in principle, this method requires one to be at the current limit of rapidity space, 
and to make stringent cuts on the transverse energy of the muons produced, which reduces 
the statistical accuracy considerably. 

An alternative approach is to use semi- inclusive DIS, in which a hadron is tagged in 
the final state in coincidence with the scattered lepton. With a deuteron target, one could 
in principle minimize the nucleon off-shell ambiguities by observing a recoil proton in the 
target fragmentation region, which would be indicative of the primary scattering having 
taken place from a near on-shell neutron [Q. Such experiments are indeed already being 
discussed in connection with Jefferson Lab and HERMES, although one must keep in mind 
the potential problem of the not being high enough at these facilities to enable a clean 
separation of the target and current regions. 

A related idea is to use semi-inclusive DIS data on hadron production in the current 
fragmentation region to measure the relative yields of 7r+ and vr" mesons. In fact, the 
HERMES Collaboration has extracted d/u from the 7r+ — tt~ difference, on both protons 
and bound neutrons [|T^. The idea here is fairly simple: at large z the u quark fragments 
primarily into a tt'^, while a d fragments into a 7r~, so that at large x and z one could have 
a direct measure of the valence d/u quark ratio. Note that one should not be too close to 
z = 1, as in that region the fragmentation process may no longer be incoherent and may 
not be factorizable into a partonic cross section and a target-independent fragmentation 
function. 

To describe the semi-inclusive process, one requires the fragmentation functions of quarks 
into pions to be known. Fortunately, the fragmentation functions for the leading m ^ vr"*" and 
non-leading u ^ 7t~ processes have been measured by the European Muon Collaboration 



at CERN W^, albeit under the assumption of charge symmetry, D'^ = , and that the 
fragmentation is target-independent. 

The apparent advantage of using both p and n targets at HERMES is that all depen- 
dence on the fragmentation functions cancels, thereby removing any uncertainty that might 
be introduced by poor knowledge of the hadronization process. On the other hand, the 
disadvantage of this approach is that one is still left with the same problem of using the 
deuteron as an effective neutron target, together with the inherent difficulties in unfolding 
the neutron structure function, which, as we argue below, are likely to be more problematic 
than the uncertainty in the fragmentation functions. The main contribution that such an 
approach could make would be to serve as a check of the extraction of the d/u ratio from 
existing inclusive DIS data on p and D. It would be preferable, given the original aim of 
avoiding nuclear uncertainties, to use only p targets, together with the empirical information 
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on the z-dependence of the fragmentation functions. 

In the parton model the number of pions in a given x and z bin can be written as 
a product of a quark distribution function, q^{x), in the hadron h, and a fragmentation 
function giving the probabihty of the scattered quark q producing a pion: 

N^-i: elq\x)Dl{z). (2) 

(The scale dependence of both the parton distribution and fragmentation functions is sup- 
pressed.) For a proton target, therefore, one has: 

Nf - 4m(x) D{z) + d{x) D{z), (3a) 
N;~ - 4m(x) D{z) + d{x) D{z), (3b) 

where D{z) = = is the leading fragmentation function (assuming isospin symme- 
try), and D{z) = D^^ = is the non-leading fragmentation function. 

In addition to protons, the HERMES Collaboration also uses deuteron targets. Using 
isospin symmetry for the parton distributions, = rf", etc., and assuming that the deuteron 
is a system of two bound nucleons, one can write for the deuteron: 

N^^ ~ (u(x) + d{x)) {W{z) + D{z)), (4a) 
Nf ~ {u{x) + d{x)) {W{z) + D{z)), (4b) 

where q denotes the smeared quark distribution of a nucleon bound in the deuteron, which 
can be approximated by: 

q(x) = [ —fN/D{y) q{x/y). (5) 

y 

The smearing function, fN/niy), gives the probability (in the infinite momentum frame) 
of finding a nucleon in the deuteron with light-cone momentum fraction y, and is given 
explicitly in Ref. [^, for example. At small x the difference between q and q is less than 
~ 2%, although at large x, where the d/u ratio is less well known, the difference can be 
much larger. 

Combining the proton and deuteron cross sections, one can define the bound neutron 
(n) distribution, = - N^: 

Nf ~ A{d{x) + e„(x)) D{z) + {u{x) + ed{x)) D{z), (6a) 
Nf ~ 4(rf(x) + e„(x)) D{z) + {u{x) + e„(x)) D{z), (6b) 

where eq{x) = q{x) — q{x). Taking the difference between the vr"*" and tt~ cross sections for 
the proton and bound neutron gives: 

Nf - Nf ~ {4u{x) - d{x)){D{z) - D{z)), (7a) 
Nf - Nf ~ (4j(x) - uix) + 4e„(x) - e,(x)) iD{z) - D{z)). (7b) 

The ratio of these. 
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P / X ^ - _ 4J(x)-72(a;) + 4e.(x)-e,(x) 

nnp[x,z) - 4u{x)-d{x) ' 

is then independent of the fragmentation function, and is a function of x only. 

Clearly, under the assumption that nuclear corrections to the quark distributions are 
neghgible, q = q, one has: 

_ 4:d{x)/u{x) - 1 

4-rf(a;)/u(a;) ' ^ ^ 

In Fig. 1 we show the ratio Rnp calculated with the assumption of no nuclear effects 
(dashed), and with the smearing correction (solid) which accounts for the Fermi motion 
and binding effects in the deuteron, as described in Ref. 0. The parton distributions were 
taken from the recent CTEQ4 parameterization |jl2| (throughout we use NLO distributions 
evaluated in the DIS scheme) at = 10 GeV^ in Fig. 1(a), while in Fig. 1(b) the d quark 
distribution was modified to have the correct perturbative QCD limit pl , according to Refs. 



d{x) d{x) 



u{x) u{x) 



+ A{x). (10) 

CTEQ4 



The analysis of Ref. suggested that at ~ 10 GeV^ the correction term could be 
parameterized by the simple form [^: 

A(x) ^ 0.2 x^ exp(-(l - xf), ^ 10 GeV^, (11) 

so that in the limit x — > 1, the modified d/u — > 1/5, consistent with the expectation from 
perturbative QCD [|^. Beyond x ~ 0.7 the difference between the corrected and uncorrected 
ratios is clearly quite dramatic. Consequently a. d/u ratio obtained from such a measurement 
in this kinematic region, without correcting for nuclear effects, would give rather misleading 
results. 

As mentioned above, one can avoid the problem of nuclear corrections altogether by 
comparing data for tt"*" and vr" production on proton targets alone. Although the statistical 
accuracy will not be as good as when both protons and deuterons are used, one can expect 
them to be sufficient to obtain reliable information on the ratio. Taking the ratio of the 7r~ 
to 7r+ proton cross sections, one finds: 

RV. = it = '^Di^)/Di^) + d{x)/u{x) 

^ ' ^ - A + d{x)/u{x)-D{z)/D{z)' ^ ^ 



The fragmentation functions D and D, measured by the EMC |TT|, are shown in Fig. 2, 



together with parameterizations for the leading and non-leading functions: 

D{z) = 0.7 {1- zy-'^^z, (13a) 
D{z) = ^^^D{z). (13b) 
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In the limit z ^ 1, the leading fragmentation function clearly dominates, D{z) ^ D{z). 
In this case the ratio {l/A)d/u [0. Although the point z = 1 cannot be reached 

experimentally, for reference we show in Fig. 3 the theoretical ratio expected in this limit 
for two different parameterizations, the standard CTEQ4 distributions [0, and with the d 
quark distribution modified according to Eq.(|10|), both at = 10 GeV^. The z = 1 limit 
provides the maximal difference between the two ratios. 

In the realistic case of smaller z, the D/D term in Eq.(|T^) will contaminate the yield 
of fast pions originating from struck primary quarks, diluting the cross section with pions 
produced from secondary fragmentation picking up qq pairs from the vacuum. Nevertheless, 
one can estimate the yields of pions using the empirical fragmentation functions in Eq.(|l^). 
Integrating the differential cross section over a range of z, as is more practical experimentally, 
the resulting ratios for cuts of 2; > 0.3 and z > 0.5 are shown in Fig. 4. One sees that 
decreasing the lower limit for z has the effect of raising the cross section ratio significantly, 
because of the larger integrated contribution from the non-leading fragmentation, which is 
more important at smaller z. Although the relative difference between the ratios for the two 
forms of asymptotic d/u behavior then becomes smaller, the absolute difference between 
these remains relatively constant, and should be measurable with the high luminosities 
available at current facilities. 

For some final words of caution, we should note that the assumptions of factorization and 
neglect of higher twists are more questionable in the very large x region. As one approaches 
the elastic limit at x = 1, the effective value of where a leading twist treatment is likely 
to be sufficient will increase with x like Q^/(l — x). From inclusive DIS data, where scaling 
sets in already at ~ 2 GeV^, one expects a leading twist analysis for = 10 GeV^ to 
be valid up to x ~ 0.8. For semi-inclusive DIS, however, it is not known up to which value 
of X the leading twist perturbative QCD framework can be applicable at any given Q^, and 
ultimately this question can only be answered experimentally. One may expect that the 
effects would not be as important for the ratios of cross sections discussed in this paper as 
for the absolute cross sections. To this effect we have also calculated at = 100 GeV^, 
and indeed find the differences with the results in Fig. 4 to be rather small. 

In summary, we reiterate the importance of an accurate experimental determination of 
the behavior of the valence d/u ratio as a; — > 1. Not only are the present fits to world data in 
clear disagreement with the predictions of perturbative QCD (unlike the reanalysis of Ref. 
0]), but the discrepancy is extremely important when it comes to estimating event rates for 
charged current events at the large values of x and accessible at HERA [|T^]. Our analysis 
of binding and Fermi motion corrections in the deuteron suggests that it is not profitable to 
try to resolve this issue by combining semi-inclusive data on protons and deuterons. On the 
other hand, the issue could be decided through semi-inclusive pion production measurements 
on the proton alone. While this will require a dedicated experiment in the kinematically less 
favored regime where the energy transfer to the final pion is larger than z ~ 0.5, in view of 
the importance of settling this question we believe it is crucial that the effort be made soon. 
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FIG. 1. Ratio Rnp in Eq.(|8|) calculated with (solid) and without (dashed) smearing corrections. 
The u and d distributions were taken (a) from the CTEQ4 parameterization and (b) with the 
d quark distribution modified as in Ref. Q to have the correct perturbative QCD limit 
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FIG. 3. Theoretical ratio, , as a function of x for fixed z = 1. The dashed line represents the 
ratio constructed from the CTEQ4 parameterization while the solid includes the modified d 
distribution according to Eq.(|lO|). 
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FIG. 4. Ratio as function of x, integrated over z between 0.3 < z < 1 and 0.5 < z < 1. The 
solid and dashed curves are as in Fig.3. 
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